Self-Doping of Gold Chains on Silicon: A New Structural Model for Si(lll)5x2-Au 
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A new structural model for the Si(lll)5x2-Au reconstruction is proposed and analyzed using first- 
principles calculations. The basic model consists of a "double honeycomb chain" decorated by Si 
adatoms. The 5x1 periodicity of the honeycomb chains is doubled by the presence of a half-occupied 
row of Si atoms that partially rebonds the chains. Additional adatoms supply electrons that dope 
the parent band structure and stabilize the period doubling; the optimal doping corresponds to 
one adatom per four 5x2 cells, in agreement with experiment. All the main features observed in 
scanning tunneling microscopy and photoemission are well reproduced. 

PACS numbers: 68.43.Bc,73.20.At,68.35.Bs,81.07.Vb 



Physical realizations of a one-dimensional metal are 
rare, in part because they may be preempted by a metal- 
insulator Peierls transition [l| . An escape clause is avail- 
able for metallic chains adsorbed on rigid substrates, 
however, since the energy penalty for the pairing dis- 
tortion may be prohibitively high. For example, when 
gold is adsorbed on silicon a variety of chain-like struc- 
tures are formed, some with unusual electronic proper- 
ties suggestive of a one-dimensional metal 0, E|- Pho- 
toemission data from the vicinal surfaces Si(553)-Au and 
Si(557)-Au reveal fractionally filled bands with strongly 
one-dimensional character, which are believed to origi- 
nate from Au "chains" just one atom wide U Structural 
models for Si(557)-Au were recently proposed based on 
total-energy calculations || and x-ray data Q. 

Even more widely studied is the parent flat surface, 
Si(lll)-Au. First reported over twenty-five years ago 
7, 8], the Si(lll)5x2-Au reconstruction has been ex- 
tensively characterized by low-energy electron diffraction 
[10, scanning tunneling microscopy (STM) [Tol HI . 
x-ray diffraction 0] , reflectance anisotropy spectroscopy 
[lij . angle-resolved photoemission 0,0], inverse pho- 
toemission , and core-level spectroscopy . Despite 
this scrutiny, its structure remains unknown. The data 
provide several constraints on any model of Si(lll)5x2- 
Au. (1) STM shows the surface to be decorated by bright 
protrusions with apparent height ~1.5 A and whose cov- 
erage, although variable, has a preferred value of one 
per 5x8 supercell 0. (2) Away from the protrusions, 
STM images show a "Y" -shaped feature whose orienta- 
tion is determined by the underlying lattice [III H2I . 
(3) Photoemission finds a strong surface band beginning 
at the the 5x2 zone boundary and dispersing downward 
toward the 5x1 zone boundary [lfij |. (4) The nature of 
this band changes from one-dimensional at the top of the 
band to two-dimensional at its bottom [i(J ■ 

In this Letter a new model is proposed for Si(lll)5x2- 
Au that explains all of these observed features. The 
model is related to those proposed earlier for Si(557)- 
Au, suggesting that all of the reconstructions formed by 
adsorption of Au on Si form a family. It is also closely 



related to the "honeycomb chain-channel" (HCC) model 
now widely accepted as the structure of the adsorbate- 
induced induced Si(lll)3xl and Si(lll)3x2 reconstruc- 
tions 21], and thus helps to unify a wide class of 
adsorbate-induced reconstructions based on Si(lll) and 
its vicinals. The present model achieves its stability 
through an unusual "self-doping" mechanism that may 
be relevant to other Au-induced Si reconstructions. 

The model is shown in Fig. 1. The reconstruction oc- 




FIG. 1: (color) Proposed "double honeycomb chain" struc- 
ture of Si(lll)5x2-Au. Large circles are Au, small circles are 
Si. The elementary 5x2 unit cell is outlined. Each unit cell 
contains two honeycomb chains (HC) based on the outlined 
hexagons, one of alternating Au and Si atoms, the other of 
all Si. Three additional Si adatoms, with 5x4 periodicity, are 
also shown (see discussion). 
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curs purely in the surface layer, and has the basic struc- 
ture of a "double honeycomb chain" (DHC) with under- 
lying 5x1 periodicity One chain is formed by hexagons 
of alternating Au and Si atoms, and one by hexagons of 
all Si (as in the HCC model). The outer Si atoms of the 
Au-Si chain are too far from the Si atoms in the Si chain 
to bond directly. As a result, two variants of the basic 
model are plausible, (i) The insertion of an additional 
"rebonding" row of Si atoms can bridge the gap between 
the chains, but only at the cost of overcoordinating the 
Si atoms in the Au-Si chain. This will be called the "5x1 
variant." (ii) Removing every other of these rebonding Si 
atoms (as shown in Fig. 1) relieves the overcoordination 
but now leaves dangling bonds in the Si chain; this will 
be called the "5x2 variant." These two variants are very 
close energetically. It will be shown below that the pres- 
ence of additional Si, occurring as adatoms, acts to dope 
the 5x2 variant with electrons and thereby to reduce its 
surface energy relative to the 5x1 variant. The optimal 
doping level occurs for one adatom per four 5x2 cells, 
identical to the observed equilibrium adatom coverage. 

First-principles total-energy calculations were used to 
determine the equilibrium geometries and relative sur- 
face energies of the basic model and its variants. The 
calculations were performed in a slab geometry with up 
to six layers of Si plus the reconstructed surface layer 
and a vacuum region of 8 A. All atomic positions were 
relaxed until the total energy changed by less than 1 meV 
per 5x2 cell; the bottom Si layer was passivated by hy- 
drogen and held fixed. Total energies and forces were 
calculated within the generalized-gradient approximation 
to density-functional theory using projector-augmented- 
wave potentials, as implemented in VASP [22I Ejj]. The 
plane-wave cutoff (180 eV) and sampling (2x4) of the 
surface Brillouin zone were sufficient to converge relative 
surface energies to within 1 meV/A 2 , adequate for the 
comparisons presented below. All the models considered 
here have equal Au coverage, and hence the relative sur- 
face energies (calculated as in Ref. EH) do not require a 
choice of Au chemical potential. STM images were sim- 
ulated using the method of Tersoff and Hamann [24| . 

Since its discovery, many structural models have been 
proposed for Si(lll)5x2-Au. Most of the early ones 
are not compatible with newer STM data, but two later 
models — from Marks and P lass (MP) [25T ] and Hasegawa, 
Hosaka, and Hosoki (HHH) [26j — have recently been 
studied theoretically [2jj ■ Neither was found to be con- 
sistent with STM or ARPES data, despite being locally 
stable with nearly equal surface energies (to within 0.1 
meV/A 2 ). The 5x2 DHC model proposed here is more 
stable than the MP and HHH models by 20 meV/A 2 , or 
2.6 eV per 5x2 cell. This energy difference is sufficiently 
large to rule out the MP and HHH models on energetic 
grounds alone. 

Within the family of the two DHC variants described 
above (5x1 versus 5x2, each with different Si adatom 



coverages) the relative surface energy changes are very 
much smaller, of order 1 me V/A 2 . This is near the lim- 
iting precision for DFT surface energies, suggesting the 
need for a simpler model to describe, for example, the 
variation of surface energy with Si adatom coverage. In 
the following, this model is developed by analyzing within 
DFT the role of the adatoms. This requires demonstrat- 
ing, first, that the observed bright protrusions in STM 
are, in fact, Si adatoms and not something else; and sec- 
ond, that the adatoms indeed supply electrons, which 
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FIG. 2: (color) Simulated filled-state STM images for double 
honeycomb chain model, (a) Image for 5x2 model (sample 
bias —0.8 V), showing "Y"-shaped feature observed in Ref.llll 
(b) Same surface, with adsorbed Si atoms as in Fig. 1 (sample 
bias —2.0 V), showing bright protrusions observed in Ref . Hfll 
Inset: simulated image from an adsorbed Au atom at the same 
in-plane location (marked by arrow), (c) Linescan through 
two bright protrusions, as marked in panel (b). 



then dope the parent band structure. 

Away from the protrusions, filled-state STM images 
show a series of side-by-side "Y" -shaped features with 
5x2 periodicity [U [H, Ei A simulated STM image 
for the 5x2 DHC model with no adatoms is shown in 
Fig. 2(a). A similar feature "Y" -shaped feature is found 
with paired "arms," a single "tail," and the same crys- 
tallographic orientation as found experimentally. 

To identify the bright protrusions seen in STM images, 
both Si and Au adatoms were considered as potential 
candidates. Binding energies for individual adatoms were 
calculated for eight possible sites on the undecorated 5x2 
surface. For Si the most favorable site is at the center of 
a Au-Si hexagon, as shown in Fig. 1, where the binding 
energy is 4.6 eV. This is much larger than the next best 
site (by 0.7 eV) and hence rules out other possible lo- 
cations. The simulated STM image, shown in Fig. 2(b), 
from a 5x4 arrangement of Si adatoms decorating these 
sites is in excellent agreement with the atomically re- 
solved images of Ref. [l9| . In particular, the bright spots 
are correctly positioned in the middle of the underlying 
row structures. The linescan in Fig. 2(c) shows their ap- 
parent height to be 1.5 A, in good agreement with the 
results of Ref. ^3- Finally, Au adatoms can be easily 
ruled out as plausible candidates: although their binding 
energies are substantial (from 3.1 to 3.9 eV), they relax 
well into the surface layer and thus produce no detectable 
STM spot, as demonstrated in the inset to Fig. 2(b). 

Photoemission data for Si(lll)5x2-Au provide a very 
stringent test for any structural model. For wavevec- 
tors along the chain direction, the ARPES data reveal a 
strong surface band beginning at the 5x2 zone boundary 
(the A2 point), and dispersing downward to its minimum 
at the 5x1 zone boundary (the Ai point) before turning 
back up 0|- The calculated electronic structure, shown 
in Fig. 3, reveals just such a band between A 2 and Ai, 
whose width and effective mass (0.6 eV and 0.4m e ) are 
in reasonably good agreement with the data (0.9 eV and 
0.5m e ). Additional bands with less pronounced surface 
character are also found in the data, and can be tenta- 
tively identified with calculated bands marked in Fig. 3. 

One remarkable feature of the ARPES results for the 
strong surface band is the continuous transition, within 
a single band, from a one-dimensional state (at the band 
maximum) to a two-dimensional state (at its minimum) 
[20| . This feature is well reproduced in the calculated 
band structure, expanded views of which are shown in 
the two small panels of Fig. 3 for wavevectors perpendic- 
ular to the chain direction. At the maximum of the strong 
band, near A2, the perpendicular dispersion is small (0.04 
eV), while at its minimum, near Ai, the perpendicular 
dispersion is much larger (0.10 eV). The ARPES mea- 
surements give very similar results, 0.03±0.03 eV and 
0.14±0.03 eV, respectively [13 . 

It is clear from Fig. 3 that the calculated electronic 
structure of the bare Si(lll)5x2-Au reconstruction is 




FIG. 3: (color) Band structure of 5x2 double honeycomb 
chain model with no adatoms. The size of each circle reflects 
the surface character of the state (radii are proportional to 
the total charge in spheres around surface-layer atoms). The 
shaded states form bands that are detected in photoemission 
data (see discussion). Left panel and right panels show disper- 
sion along and perpendicular to chain direction, respectively. 
The colored circles in the left and right panels mark the same 
states; note the change of energy scale. 

metallic. As suggested in Ref. Q, it is interesting to 
ask whether the addition of Si adatoms could render 
the system either insulating, or at least "less metallic"; 
whether it is energetically favorable to do so; and — if it 
is — whether the predicted optimal coverage of adatoms 
corresponds to the equilibrium coverage deduced from 
STM data, one per 5x8 supercell ^3- The answer to all 
three questions is affirmative, as shown below. 

Near the Fermi level the bands are very close to one- 
dimensional, as discussed above. Hence, the number of 
additional electrons required to render the band struc- 
ture (provisionally assumed to be rigid) insulating can be 
trivially determined from Fig. 3 to be two per 5x2 cell. 
To determine what coverage of Si adatoms is required to 
provide two electrons per 5x2 cell, explicit band struc- 
ture calculations must be performed. For a coverage of 
one adatom per 5x2 cell, the resulting band structure is 
insulating. This implies that the addition of one adatom 
to the bare 5x2 surface creates only one additional state 
in the occupied manifold: two of the adatom's valence 
electrons fill this state, and the other two dope the par- 
ent band structure, rendering the full system insulating. 
These conclusions are independent of any assumptions 
about the rigidity of the parent band structure. Never- 
theless, it is important to note that the bands do shift 
quite rigidly, with only small changes of order 0.1 eV or 
less. For lower adatom coverages, smaller but similarly 
rigid shifts are expected. This transforms the determina- 
tion of the optimal adatom coverage into the equivalent 
task of computing the optimal electron doping. 

In the absence of adatoms, the DFT surface ener- 
gies of the 5x1 and 5x2 DHC models are very close, 



4 




0.00 0.25 0.50 0.75 1.00 
Doping (electrons/5x2-cell) 



1.25 



FIG. 4: Variation of surface energy with electron doping, for 
5x1 and 5x2 models. The most favorable structure has 5x2 
periodicity with 0.5 extra electron per 5x2 cell. 



with the 5x1 variant preferred by less than 1 meV/A 2 . 
Changes in the surface energy due to electron doping 
were modeled by calculating the DFT total energies for 
cells with additional electronic charge (with a compen- 
sating background charge to preserve overall neutrality, 
plus the standard correction to treat the resulting spu- 
rious interactions (HI). The results are shown in Fig. 4. 
For the 5x1 variant, any additional electronic charge in- 
creases the surface energy above its undoped value. The 
5x2 variant behaves quite differently: its surface energy 
is minimized for a doping level very close to 0.5 elec- 
tron per 5x2 cell. Since each adatom was earlier shown 
to supply two doping electrons per 5x2 cell, this opti- 
mal doping level can be most easily achieved with one 
adatom per four 5x2 cells. Moreover, at this optimal 
adatom coverage the surface energy of the 5x2 variant is 
lower than that of the 5x1, suggesting that this is indeed 
the observed ground-state equilibrium phase j2t|. These 
findings are in excellent agreement with the experimental 
observation that when excess Si is evaporated at low tem- 
perature onto a Si(lll)5x2-Au surface with one adatom 
per four 5x2 cells, subsequent annealing at higher tem- 
perature will cause the extra Si to diffuse away and return 
the system to its equilibrium state. 

In summary, a new structural model has been pro- 
posed for the Si(lll)5x2-Au surface, consisting of a 
"double honeycomb chain" reconstruction decorated by 
Si adatoms. Simulated STM images from this model re- 
produce a number of experimentally observed features, 
including the in-plane location and apparent height of 
the bright protrusions due to the adatoms, and the "Y"- 
shaped features seen away from these protrusions. The 
calculated band structure reproduces the main features of 
recent photoemission data, including the unusual change 
in dimensionality observed within the main surface band 
between its energy extrema. Finally, the Si adatoms act 
as electron donors that dope the parent 5x2 band struc- 



ture, reducing its surface energy and stabilizing it relative 
to other models. The optimal doping level is equivalent 
to one Si adatom per four 5x2 cells, in agreement with 
experimental observation. 
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